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Abstract 

The scalar potential, time component vector potential and flux tube quark confinements are 
studied in this paper. We find that the predictions of scalar confinement and time component vector 
confinement are in considerable conflict with measured values while the flux-tube confinement works 
well to explain the experimental data. We also study the relativistic flux tube model. From the 
comparison of the exact numerical solution with the analytic approximation solution for heavy- 
light mesons, we find that the solutions are much more in agreement with each other for higher 
excited states since the deep radial limit is better satisfied. 
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I. INTRODUCTION 



The scalar potential quark confinement has been proposed ctS ct SI gnificant model in 
hadron physics. In order to account for the observations of spin-orbit splitting, that the 
short range Coulomb potential was a Lorentz vector and the long range confining potential 
a Lorentz scalar was proposed by Henriques, Kellett, and Moorhouse There is no color 
magnetic field to act on the quarks' spin. As Buchmiiller pointed out that a color electric 
flux tube should automatically yield the desired pure Thomas spin-orbit interaction because 
there is no color magnetic field in the quark rest frame. The Milan group clarified this 
question dramatically around 1990. Using the low velocity Wilson loop formalism pioneered 
by Eichten and Feinberg and later by Gromes 0, they found both the spin- dependent 
and spin- independent relativistic corrections in heavy oniaj^]. Their result was that the 
long-range spin-independent QCD corrections differed from those of scalar confinement, 
though the long-range spin-orbit corrections were the same pure Thomas ones given by 
scalar confinement. These results were subsequently verified by lattice simulations of QCD 
0. The scalar quark confinement models remain popular in spite of the inconsistency 
with QCD, probably because of their relative ease of solution. We point out here that the 
predictions of scalar confinement also conflict directly with experiment. We study the scalar 
potential, time component vector potential and flux tube quark confinements and relativistic 
flux tube model in this paper. We investigate the dependence of heavy-light meson mass 
on the mass of the light quark. This paper is organized as follows. We exhibit a solution 
for scalar quark confinement in section II. In section III, we compare the predictions of the 
quark confinement with the experimental result and find that the predictions of the scalar 
confinement and the time component vector confinement are in considerable conflict with 
measured values while the flux-tube confinement works well to explain the experimental 
data. We review the relativistic flux tube model in section IV. Under two limit we solve 
the equations for the relativistic flux tube model in section V. An exact numerical solution 
is also demonstrated . In section VI, the comparison of the numerical solution with an 
analytical approximation solution is discussed. The conclusions are drawn in section VII. 

II. SOLUTION OF THE SCALAR QUARK CONFINEMENT EQUATION 

We consider a potential field consisting of central Lorentz scalar and time component 
vector fields, S(r) and V(r) respectively. To simplify our derivations, We assume that 
scalar confinement S(r) = a\r\ = ar and a time component vector short range interaction, 
U(r) = — A = —-. For a spinless quark m moving in this potential field, we find 



We wish to treat the short range parameter k and the light quark mass m as small pertur- 
bations. In the case where both k and m vanish, the zeroth order Hamiltonian is 




(1) 




(2) 



where 
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The Hamiltonian Ho has the same eigenstates as its square which, with the replacement 

id 2 

r dr 2 



Pr^~ -4-2 r i leads to the harmonic oscillator equation 



cPuo ^ £(£+1) 2 2 \ _ m 

-^2 + [ E o ~2 ar Ju = 0, (4) 

where in the limit of small r, uq — > r e+l and Uq is normalized to 

dr | wo | = 1. (5) 
The solution for the wave function and eigenvalue is standard; 

Mr) = M n/ r e+1 e- l2ar2 d(ar 2 ), (6) 

^ " r(£ + n + i) ' (7j 

El = 2a (l + 2n-l). (8) 
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Here L n _\(ar 2 ) is the usual Laguerre polynomial, n is a positive integer starting with 1, and 
£ is a non-negative integer starting with zero. 

We now determine the effect of turning on k and m using the Feynman-Hellmann theorem 

dE/dX = (dH/dX), (9) 

where A is any parameter of the Hamiltonian. Taking the expectation values using the 
k = m = wavefunctions © will yield an expansion in these parameters. To leading order 
we have 

„ „ / 1 \ ma , . , 
E = E -k (r- l ) + — (r) + ... . (10) 

The expectation values are worked out in general in the Appendix. Here we only consider 
the m and k dependence of the IS and IP states (i.e., n — 1 and £ = and 1). The results 

are, 

(r-\ s = (11) 

(r-\ P = IS < U > 
(r) is = ~L> (13) 



dE ls 2 



Using Eq. (|TU|) . we obtain the m dependence near k = m = (m slope) for the energies, 

0.6515, (15) 
1) = 0.0214. (16) 



dm V37T 
d(E 1P -E ls ) 2/4 



dm ^/37T V\/l5 
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Again using Eq. (jlOjl we obtain the k dependence near k = m = (k slope) for the energies, 



1 dE ls 
yfa dk 

1 d{E 1P -E ls ) 
dk 




= -1.128, (17) 
0.376. (18) 

Next we compare our analytic values for the slopes at m = k = in Eqs. ()15j) to (j!8J) with 
our numerical method. This step is important because a general analytic solution of the 
spinless Salpeter equation is not known. Only for the remarkable case of a massless particle 
and linear scalar confinement, which is equivalent to the non-relativistic harmonic oscillator, 
can one obtain an analytic solution. In the general case, one must rely on exact numerical 
solutions. Some time ago, a variational method, the Galerkin method, was introduced 
into particle physics to solve the spinless Salpeter equation with a time component vector 
interaction jil]. This very robust method is applicable to a wide range of differential and 
integral equations. The method has been sharpened over the years by many authors 0. 
One can cope with eigenvalue equations for operators that are complicated functions of both 
momenta and coordinates, such as the scalar confinement Hamiltonian (0), by using basis 
functions that can be Fourier transformed. We have performed a careful numerical solution 
of the eigenvalue equation for the Hamiltonian (Q), for small k and m, and found the m 
and k slopes. The results are in excellent agreement with the values obtained by analytic 
calculation in Eqs. (fT3j) through (fTHj) . 



III. COMPARISON OF THE PREDICTIONS OF DIFFERENT QUARK CON- 
FINEMENTS WITH EXPERIMENTAL RESULTS 

A. Spin splitting is independent of light quark mass 

We first use experimental data to demonstrate, rather conclusively, that spin splittings 
within a given orbital angular momentum multiplet do not depend on the light quark mass. 
In particular, we consider several heavy-light meson spin multiplets in which both the strange 
and non-strange members have been observed [ljj. First we examine the D and D s type 
mesons with i = (<?- waves). The hyperfine splittings for D s and D mesons are 

D* — D s = 143.8 ± 0.4 MeV, 

D* ± -D± = 140.64 ± 0.10 MeV, (19) 
D* - D = 142.12 ± 0.07 MeV. 

The corresponding S'-wave hyperfine splittings for B type states are 

B* - B s = 47.0 ± 2.6 MeV, 

B* - B = 45.78 ± 0.35 MeV. (20) 

It is clear that the substitution of a strange quark for a non-strange light quark changes the 
.S'-wave hyperfine differences by at most a few MeV. 
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We next consider some measured P-wave heavy-light spin splittings. From 10] we find, 

D s2 -D sl = 37.0 ±1.6 MeV, 
Dl - P>° = 36.7 ± 2.7 MeV, (21) 
Df - Df = 32 ± 6 MeV. 

Again we note the apparent vanishing of light quark mass dependence this time in a 
P-wave spin splitting. We conclude that spin splittings are only weakly dependent on the 
light quark mass. 



B. The light quark mass dependence of a IP-IS difference 

In the preceding subsection we observed from experiment that both S'-wave and P-wave 
heavy-light spin splittings (within a spin multiplet) were independent of the light quark 
mass. We next consider the mass splittings between pairs of states corresponding to different 
orbital angular momenta and examine the light quark mass dependence of this difference. 
We choose the D\ P-wave state and the pseudoscalar D meson. The best measurements are 

A u = Dl - D° = 557.5 ± 2 MeV. (22) 

When the u light quark is replaced by a strange quark, the corresponding difference becomes 

A s = D sl - D s = 567.3 ± 0.4 MeV. (23) 

The differences A u and A s are amazingly similar. We see that they differ by 

A = A s - A„ = 9.8 ± 2 MeV. (24) 

We demonstrated in Sec. IIII Al that both S- and P-wave spin splittings are, within error 
and isospin uncertainty, independent of light quark mass. We may therefore conclude that 
the m-dependence of the difference A s — A u also represents the m-dependence of the spin- 
averaged IP-IS excitation energies. To show this explicitly, we write the HL meson mass 
as the sum of the heavy quark mass and the excitation energy, 

M = m Q + E. (25) 

We then separate the excitation energy into spin-averaged and spin- dependent parts, 

E = E SA + E SD . (26) 

As in Eqs. (I22J) and (}2*3*|) . we define A s and A u to be the differences between P-wave (J p = 
1 + ) and S*-wave (J p = 0~) states, for strange and non-strange light quarks respectively. 
Using our observation that the spin-dependent parts are essentially independent of light 
quark mass, we see that the spin dependent parts will cancel in the difference A = A s — A u . 
Thus 

A = (E 1P - E 1S ) S - (E 1P - E 1S ) U (27) 

measures the light quark mass dependence of the spin-independent IP-IS" level difference. 

Two main results emerge from this subsection. First, the light quark mass dependence 
is only about 2% of the measured level difference. This is a useful tool in the analysis of 



heavy-light spectroscopy [11[. Second, although it is small, A is definitely not zero. We will 
use the actual value (|2^j) to test the predictions of various assumptions about the nature of 
quark confinement. 
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C. Comparison with experimental data 

We are now prepared to compare carefully the different confinement predictions with 
experiment. As we noted, the m dependences of all the heavy-light states are amazingly 
similar. We have also noted that this universal m dependence is nearly satisfied in an 
analytical calculation. In the example with scalar confinement, the m slope for the difference 
IP— IS is about 30 times smaller than each separate slope. Furthermore, we noted in 
Sec. IIII Bl that when one compares two heavy-light states, one P-wave and one S-wave, the 
difference changes by less than 2% when a non-strange light quark is replaced by a strange 
one. This change is not zero however, but for the D\ and D states has the experimental 
value J2U) 

A = A s - A„ = 9.8 ± 2 MeV. (28) 

In FIG.l we show the numerical mass splittings of heavy-light mesons for three different 
confinement scenarios, all with the same short range energy —k/r, and the same k = 0.5, 
as a function of light quark mass. Each of the three confinement scenarios has the same 
asymptotic confinement force, a = 0.18 GeV 2 . The upper curve assumes linear scalar 
confinement, the middle curve is the prediction of the relativistic flux tube and the lower 
curve linear time component vector confinement. In the scalar and time component vector 
potentials, the potentials S(r) and the long-range part of V(r) respectively are ar. In the 
case relativistic flux tube model the string tension is a. 

We note that the scalar confinement has the most rapid increase of the IP— IS difference 
as the light quark mass increases. The flux tube confinement is intermediate and time 
component vector linear confinement actually decreases slightly. Using the reasonable values 
for the strange quark mass (500 MeV) and the non-strange mass (300 MeV), we find the 
following values for A, 

A sca i ar = 19 MeV, (29) 

^flux tube 

10 MeV, (30) 
A TC v = -1 MeV. (31) 

In comparing these values to the experimental value given in Eq. (1231) . 

A cxp = 9.8 ± 2 MeV, (32) 

we observe that both scalar confinement and time component vector confinement are quite 
inconsistent with the experimental result. However, the flux tube confinement prediction is 
in agreement with the experimental result. 

IV. REVIEW OF RELATIVISTIC FLUX TUBE EQUATIONS 

From the above study on the different quark confinement models, we observe that flux 
tube confinement prediction is in agreement with the experimental result. Next we continue 
our research on the flux tube model and concentrate on the relativistic flux tube (RFT) 
model|I3. The RFT model or QCD string play an important role in describing the rela- 
tivistic meson states. Its fundamental assumption is that the QCD dynamical ground state 
for large quark separation consists of a rigid straight tube-like color flux structure connecting 
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FIG. 1: Dependence of the IP— IS energy level difference on the light quark mass m. The numerical 
results compare different confinement mechanisms (all with the same long range confinement force 
a = 0.18 GeV 2 and all three calculations assume the same short range Coulombic constant k = 0.5.) 



the quarks and for massless quarks it reduces to Nambu string. jiff] This model is consis- 
tent with both spin-dependent and spin-independent QCD expectations while the scalar 
confinement potential model0] is not consistent with spin-independent QCD relativistic 
correction. Even though with this simple physical picture the RFT model coincides with 
the QCD relativistic correction derived from the rigorous Wilson loop formalism. We now 
review the spinless quarks and straight QCD string model and give the RFT equations. 
Consider a heavy-light (HL) meson made up of two spinless quarks, one heavy quark with 
mass M and one light quark with mass mi. They satisfy the conditions as M > mi -> 
and M ^> ar where a is the confinement constant and r is the interquark separation. Since 
the center of mass remains close to the heavy quark in the HL meson, the constraint that 
the total momenta perpendidicular to the rotation axis and to the interquark axis vanish 
is automatically satisfied. For a heavy-light meson (HL), we obtain the equations of QCD 
string as [Hi]. 



h. = W rl± v ± + 2arf(v ± ) (33) 
E = W r7± + arS(v ± ) + V(r) (34) 
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whereL and E are angular momentum and energy of the meson respectively, V(r) repre- 
sents the short range potential. We also define the functions: 



S(v ± ) = aTCS ^ (35) 

1 



fM = —[SM- y/l-vl] (36) 



W r = \jv 2 r + m 2 (37) 
7l = T^-T (38) 



v 



_L 



V. SOLUTIONS UNDER TWO LIMITS 



A. Circular Motion limit 



We now consider the light quark of a HL meson moves circularly around the heavy quark 
which is at rest in the center of the circular motion. So this will reduce the p r — > and 
v± — > 1. This reduction yields 

S(v ± ) = \ (39) 
fM = I (40) 

When pluging into the general HL RFT equations, we obtain 



^ = ^ (41) 
r 4 

„ nar , , 

S = — (42) 



Therefore we find the l Nambu' slope for the HL meson is 



L 1 

«™ = = ™ (43) 



B. Deep radial limit 



That the QCD string reduces to a time component vector potential was shown under the 
deep radial limit in which L E 2 . It is reasonable to apply the semi-classical quantization 
condition for a spherically symmetric system to quantize the QCD string, 

1 r r + 

(n + -)tt = / drp r = I (44) 

2 Jr_ 
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where r+, r_ represent the radial distances at the turning points of the motion. Using 
this semi-classical quantization condition we obtain the following equation satisfied by the 
semi-classical TCV integral for HL meson, 



where 



2(e 2 -L)J(e 2 + L) l 

^ ) -J=(n + ±)« (45) 



3 = I dyl — i^w — (46) 



m = h = — e 2 = — 47 

e 2 + L e 2 4a 



We have defined that 



When evaluating Eq. ( I45j) we need replace the classical angular momentum with Langer 
correction to take into account the centrifugal singularity, namely L — > I + 1. 

We give an exact result of this TCV integral for the heavy-light mesons by means of 
numerical method. The relation(|43j) of angular momentum and energy in ground state of 
the system when k = is shown in FIG. 2. 

To comprehend more about the variations of angular momentum with energy in different 
state of the quark system, we draw three states in FIG. 3 corrresponding to n = 0, 1 and 2 
respectively, 

Now we are ready to turn on the color Coulomb potential V(r) = — - to explored what 
effect the short range Coulomb potential has on our graphs. We find the new relation of 
angular momentum and energy in the system is 



f ^ 2 (e 2 -L + k)^(e 2 + L + k) 

(n + -)tt = Jt (48) 




According to the above equation, we may draw the graphs representing the relation of 
angular momentum and energy numerically. Considering the ground state of the quark 
system and choosing k — 0, 0.2, 0.4 separately, the modified curves due to the short range 
color Coulomb potential are shown in FIG. 4. Note that for a fixed L, higher k value curve 
gives the lower energy. This is due to the color Coulomb potential we take. 



VI. COMPARISON OF AN ANALYTIC APPROXIMATION SOLUTION WITH 
THE EXACT NUMERICAL SOLUTION 

A. An analytic approximation solution 

We now derive an analytic approximation of RFT considering L E 2 and k = 0. To 
this purpose, we define 
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FIG. 2: Dependence of L on e 2 in the ground state of a heavy-light meson. 



Since (3 <C 1 in the deep radial limit, m and h may be represented with f3 as 



in = r ^ ~ 1 — 2/3 



1 + 



/i« 1-/3 



Plugging into Eq.([46|) we obtain 



Assuming that, 



-1 ) /(l- y 2)(l-y2 + 2/?y2) 

J =' " y T^TT^? 



1 +X 



then J may be represented with the new parameters x and /3 as 



10 



L 



3 

2.5 

2 

1.5 

1 

0.5 



/ 



2 3 4 5 6 

FIG. 3: The middle dashed line is for the first excited state, the right one is for the 2nd excited 
state while the solid line is for the ground state. 



2 Jo 1 + 



1 + 2(3x) 



x 



(54) 



io i + (3x 

Next expand J with taylor expansion and drop terms of second order in (3 or higher, so 



For = 0, 



and 



<9 T 



1 r°° 
J(0) = - / dx 



2 Jo \x{l + x)- : 



(55) 



(56) 



dJ_ 

d(3 



1 roo 



f3=0 



2 Jo 



1 r c 
' 2 Jo 



dx\ 



dz- 



y/l + 2px s 
1 + px ' 



x(l + x) 3 

1 



0(3 



7T 



v / aT^)(i+,) 2 2 4 

where we let z = (3x. Therefore when (3 1, we find to first order in (3 



(57) 
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FIG. 4: The solid curve represents color Coulomb constant k = 0, the left dashed curve corresponds 
to the k = 0.4 while the middle dashed one is for k = 0.2. 



/ = 2(62 " L)V?TI J = 2e'(l - + + (i - 1)0] 



2 4 



7T 



= 2^(1 - -/?) 



(58) 



After substituting e 2 with j-, we obtain an analytical approximation of the heavy-light 
semi-classical time component vector potential integral, 

Applying the semi-classical quantization condition for a spherically symmetric system 
and replacing the classical angular momentum with Langer correction to take into account 
the centrifugal singularity, we immediately obtain the spectroscopic relation 



L + 2n + 



3 _ & 
2 ixa 



(60) 
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FIG. 5: The ground state n=0 
which is an excellent approximation to its numerical solution shown in section IV 1BI 



B. Consistency of the analytic approximation solution and the exact numerical 
result 

We are now ready to compare this analytic approximation result with the exact numerical 
one. We draw the six figures(from FIG. 5 to FIG. 10) for every other n starting from to 10. 
From these figures we easily identify that the approximations are more accurate along with 
the increasement of n since the deep radial limit is much more satisfied. 

VII. SUMMARY AND CONCLUSION 

We have studied the scalar potential, time component vector potential and flux tube 
quark confinements. The Lorentz scalar quark confinement has a long history and is still 
widely used despite its well-known theoretical faults. We point out here that the predictions 
of scalar confinement also conflict directly with experiment. We investigate the dependence 
of heavy-light meson mass differences on the mass of the light quark. In particular, we 
examine the strange and non-strange D mesons. Under two limit we also solve the relativistic 
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FIG. 6: The 2nd excited state 

flux tube equations. We have presented the exact numerical solution for HL mesons and 
compared the numerical solution with an analytical approximation solution. 

The conclusions are that spin splittings are independent of the light quark mass value. 
We then can extract the change (|24|) in the spin-averaged IP— IS level difference as one 
replaces the non-strange by a strange light quark. We compare the predictions of the scalar 
potential, time component vector potential and flux-tube quark confinement scenarios with 
experimental results. We conclude that flux-tube confinement works well while both scalar 
and time component vector confinement fail badly to explain the experimental data. We 
observe therefore that scalar confinement has at least one specific point of disagreement with 
experiment. This complements the theoretical disagreements with QCD mentioned in the 
introduction. For the relativistic flux tube model, from the comparison of exact numerical 
solution with the analytic approximation solution for heavy-light mesons, we see that they 
are more in agreement with each other for higher excited states since the deep radial limit 
is better satisfied. The relativistic flux tube model gives the correct universal relations for 
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FIG. 7: The 4th excited state 



HL meson spectroscopies. 



Acknowledgment 

The author would like to thank M. Olsson, C. Goebel and S. Veseli for their helpful 
discussions and collaborations. This work was supported in part by the U.S. Department of 
Energy under Contract No. DE-FG02-95ER40896. 



15 



L 



. 9 



. 8 



. 7 



. 6 



11.1 11.2 11.3 

FIG. 8: The 6th excited state 
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Appendix Derivation of the expectation values 

We now derive the expectation value of r p with the harmonic oscillator wavefunctions 



(r p ) = Nl e / drr 
' Jo 



2£+2+p e ~ar 2 



A change of integration variable to the dimensionless combination z = ar 2 yields 



3 i P 



2a" + ^(r p ) =NZ, / dzz l+ ^e 
' Jo 



(61) 



(62) 
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FIG. 9: The 8th excited state 



It is helpful to use the Chu-Vandermonde sum formula j!5| . 
where the Pochhammer symbol (z)n is defined as, 



(z) N = z{z + l)---(z + N-l) 

Wo 



T(z + N) 
T{z) 
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FIG. 10: The 10th excited state 



With the choices 

(66) 
(67) 

n - 1 - ^. (68) 

we substitute Eq. (J63|) into (J62)) and use the orthonormality relation for Laguerre polynomials 



a = £+-, 
2' 
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dz z? e- U$(z) 4(z) = T{j + & + 1} 6 jf , 



to obtain our general result, 

(r p ) = £ 



P 

2/ n-j 



j'-ly r(n + £+|)(n-j')!' 
In Eq. (170)1 we use the notation for the binomial coefficients, 



77 



77! 



.777 / m!(77 — 777,)! 

The specific result that is required in Section 2 is for the ground state (77 = 1) is 



(r p ) . 

< ' n=l 



_£ 

a 2 



r(£ + §^ 



2 1 2- 



Some cases of direct interest are 



p = —1 
p = 0: 
p = 1 : 

J7 = 2 : 



r(i + 1) 
r(£ + 



3^ ' 



(1> = 1, 



(0 



(r 2 ) 



1 r(i + 2) 
7^r(£ + |)' 



+ 



(69) 



(70) 



(71) 
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